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In this paper we study the breaking of chiral symmetry with non-confining power-like potentials. 
The region of allowed exponents is identified and, after the previous study of confining (positive 
exponent) potentials, we now specialize in shorter range non-confining potentials, with a negative 
exponent. These non-confining potentials are close to the Coulomb potential, and they are also 
relevant as corrections to the linear confinement, and as models for the quark potential at the 
deconfinement transition. The mass-gap equation is constructed and solved, and the quarks mass, 
the chiral angle and the quark energy are calculated analytically with a exponent expansion in the 
neighbourhood of the Coulomb potential. It is demonstrated that chiral symmetry breaking occurs, 
but only the chiral invariant false vacuum and a second non-trivial vacuum exist. Moreover chiral 
symmetry breaking is led by the UV part of the potential, with no IR enhancement of the quark 
mass. Thus the breaking of chiral symmetry driven by non-confining potentials differs from the one 
lead by confining potentials. 



I. INTRODUCTION 

The problem of the spontaneous chiral symmetry 
breaking (xSB) is one of the QCD cornerstones. While 
xSB driven by a confining potential has been studied in 
detail, we now explore the effect of sorter range, non- 
confining, power-law potentials in xSB. 

Continuing the well defined mathematical problem 
of studying xSB driven by power-law potentials, [l|, 
V(r) = ±K 1+l3 r 13 , we now specialize in negative ex- 
ponents (3 < 0. For diferent values of /?, we find numer- 
ically chirally noninvariant possible vacua of te theory, 
solutions to the corresponding mass-gap equation. We 
exploit the potential model for QCD [lj, whose origins 
can be traced back to QCD in the truncated Coulomb 
gauge and which is proven to be successful in studies 
of the low-energy phenomena in QCD Q. This class 
of models can be indicated as Nambu and Jona-Lasinio 
type models Q with the current-current quark interac- 
tion and the corresponding form factor coming from the 
bilocal gluonic correlator. A standard approximation in 
such type of models is to neglect the retardation and to 
approximate the gluonic correlator by an instantaneous 
potential of a certain form. 

There are two different motivations to study xSB 
driven by a class of non-confining potentials, the correc- 
tions to xSB due to non-confining potentials, and x^B 
at the deconfinement transition. Notice that the quark- 
antiquark static potential computed in lattice QCD has 
clearly two distinct components, confinement (linear-like 
) and the shorter range Coulomb-like potential. In par- 
ticular the shorter Coulomb-like range part of the quark- 
antiquark static potential may be more complicated than 
a pure — Coulomb potential. There are at least two 
different Coulomb potentials, the pcrturbativc Coulomb 
which includes logarithmic corrections and the Luscher 
Coulomb due to the confining string fluctuations [I]- 
Moreover the matching of these two Coulomb potentials 
and of the long range linear potential may also be de- 
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FIG. 1: We show examples of non confining potentials, in 
particular the T > T c Lattice QCD data for the free energy 
Fi, thanks to 0,@,!,|I3,[Il|]Olaf Kaczmarek et al. The solid 
line represents the T — static quark-antiquark potential. 



scribed by a non-confining potential. Thus potentials in 
the neighbourhood of a pure Coulomb potential, with 
j3 7^ 1 may also be phenomcnologically relevant. More- 
over, at the deconfinement, say at the deconfinement 
phase transition of QCD or at a large number of 
flavours as in walking technicolour Q , the confining po- 
tential vanishes and it is then relevant to study the im- 
pact of potentials, shorter range than the confining po- 
tential, in the spontaneous x^B. To illustrate that dif- 
ferent potentials may be relevant, in Fig. [T] we show 
different finite temperature T free energies computed in 
lattice QCD @, H, H G3 EH by Kakzmarek et al. Notice 
that here we only address chiral symmetry breaking at 
zero T, nevertheless the present work may also be used 
as a starting point for the study of x^B a finite T or at 
finite /i, where confinement is lost. 

The problem of instability of the chirally invariant vac- 
uum for power-like confining potentials, i. c. for positive 
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exponents, has already been studied thoroughly in the 
middle of 80's by the Orsay group [H,[l3|,E3,[l5| and such 
an instability was proved for the range < (3 < 3. For 
numerical studies, the harmonic oscillator type potential, 
j3 = 2, was chosen by these authors, as well as by the Lis- 
bon group [H, [13, EH and the Dubna [l9[ group, and a 
set of results for the hadronic properties were obtained 
in the framework of the given model. Adler and Davis, 
the Lisbon group, the Zagreb group, and the Rayleigh 
gro up also studied the linear potential many years ago 

M ei, ism. 

Recently Bicudo and Nefediev [l[ solved the mass-gap 
equations for < (3 < 2 explicitly, and demonstrated that 
the chiral angle, the vacuum energy density, and the chi- 
ral condensate are smooth slow functions of the form of 
the confining potential, so that the results obtained for 
the potential of a given form - the linear confinement be- 
ing the most justified and phenomenologically successful 
choice [3, HjJ HH, H3, have a universal nature for any 
quark-quark kernels of such a type. Following the set 
of recent publications devoted to possible multiple solu- 
tions for the chirally noninvariant vacuum in QCD j23 | 
(see also [25| where a similar conclusion was made in a 
different approach), Bicudo and Nefediev also addressed 
the question of replicas existence for various power laws 
r@ , and found that for the whole range of allowed powers, 
< (3 < 2, replica solutions do exist similarly to the case 
of (3 = 2 studied in detail in [ijj ■ 

This prompted us to extend the Coulomb potential 
with other negative exponents, and to study in detail it's 
contribution to xSB. In Section II we extend the mass 
gap equation for power-law potentials with negative ex- 
ponents. In Section III we study analytically the mass 
gap equation. In Section IV we solve algcbraicly the mass 
gap equation in the chiral limit. In Section V we address 
the quark energy and the vacuum energy. In Section VI 
we conclude. 



II. THE MASS GAP EQUATION FOR THE 
POWER-LAW POTENTIALS 

We now derive the mass gap equation for the power- 
law potentials. This extends the derivation of Bicudo and 
Nefediev [l[ for positive exponents (3. 

The chiral model which we use for our studies is given 
by the Hamiltonian with the current-current interaction 
parametrized by the bilocal correlator K^,, 

H = J d 3 x^( X ,t)(-if ■ v)^(x,i) + 

\ |d 3 xd 3 y^(x,t)^(x-y)^(y,i), (1) 

where the quark current is J°(x, t) = ?/>(x, <)7 M ^-'0(x, t) 
and the gluonic correlator is approximated by a density- 
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FIG. 2: We show the dynamical masses m(p), solutions of the 
mass gap equation with zero bare mass mo = 0, obtained by 
Bicudo and Nefediev The different solutions correspond 
to the positive exponents (from left and bottom to right and 
top) (3 = 0.1, 0.3, 0.5, 0.7, 0.9, 1.1, 1,3, 1.5, 1.7 and 1.9. Notice 
that the dynamical mass vanishes in the limit of f3 — > 0. The 
masses are a function of the momentum p and dimensionless 
units of Ko = 1 are used. Also notice that the confining 
potentials enhance the masses in the IR. 



density potential, 

K$(x - y) = s^o^odx - y|), (2) 

3 

where the denominator | normalizes the Gcll-Mann ma- 
trix contribution to the mass gap equation. We now 
study the class of potentials with 

V (\r\) = -aK P +1 \rf , (3) 

where the only dimensional parameter of the model is the 
strength of the confining force Kq. 

Previously Bicudo and Nefediev [l| studied the power- 
law confining potentials, using the notation ii"o a+1 |x| a , 
adequate for the study of positive exponents, in partic- 
ular for the class of confining potentials < a < 2, in- 
cluding the linear and the harmonic oscillator potentials. 
The results of Bicudo and Nefediev are show in Fig. [5] 

Here we specialize in shorter range potentials, with a 
negative exponent. Since the Coulomb potential is fre- 
quently noted —a/r, with dimensionless a, we now adopt 
the notation of eq. ^ , where the exponent is denoted (3. 
Another difference to the previous work of Bicudo and 
Nefediev [l[ is the sign of the potential. While the con- 
fining potentials are attractive, it is necessary to have a 
negative sign for the non-confining potential for the po- 
tential to be attractive, and for the existence of bound- 
states in the spectrum. 

The relativistic invariant Dirac-Feynman propagators 
[l5| . can be decomposed in the quark and antiquark 
Bethe-Goldstone propagators [l||], used in the formalism 
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of non-rclativistic quark models, 

i 
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= -icr 2 75<(k) , 

where it is convenient to define, 

S(k) = sin^(fc) = m{k)D{k) 
C(k) = cosip(k) = kD(k) 

D(k) = 



(4) 



i 



\Jk 2 +m{kf 



(5) 



where m(k) is the constituent quark mass and (p is the 
chiral angle. In the non condensed vacuum, ip is equal 
to arctan In the physical vacuum, the constituent 

quark mass m(fc), or the chiral angle ip(k) = arctan m ^ , 
is a variational function which is determined by the mass 
gap equation. We illustrate here examples of solutions, 
for the positive exponents /? > depicted in Fig. [2J 

There are three equivalent methods to derive the mass 
gap equation for the true and stable vacuum, where con- 
stituent quarks acquire the constituent mass [26]. One 
method consists in assuming a quark- antiquark 3 Pq con- 
densed vacuum, and in minimizing the vacuum energy 
density. A second method consists in rotating the quark 
and antiquark fields with a Bogoliubov-Valatin canonical 
transformation to diagonalize the terms in the hamilto- 
nian with two quark or antiquark second quantized fields. 
A third method consists in solving the Schwinger-Dyson 
equations for the propagators. Any of these methods lead 
to the same mass gap equation and quark dispersion re- 
lation. Here we replace the propagator of cq. ([J| in the 
Schwinger-Dyson equation, 

dw' d 3 k' 
~2T (2tt) 3 ' 

u(k') s ,u^k') s , v(k') s ^(k') 6 



= u\(k) \ kk ■ a + mo/3 

E 



w'-E(k')+ie -w'-E(k')+ie 
< - r dw' d 3 k' 

2tt (2tt) 3 
tt(*').'« + (fc').' v(k') s ^(k') s , 



E(k) = ul(k) \ kk ■ a + m /3 

E 



w'-E(k')+ie -w'-E(k')+ie 




the mass gap equation and the quark energy are finally, 



= +S( P ) B(p) - C( P ) A(p) 
E(p) - +S(p)A(p) + C( P )B(p) 



(7) 

(8) 



where 



A{p) 
B(p) 



d 3 k 



V(p-k)S(k) 



(2tt) 3 
1 f d 3 k ~ 
P+ 2j (2^(P- k )(P-W) (9) 



Using the Bogoliubov-Valatin transformation, the Hamil- 
tonian (TTJ) splits into the vacuum energy, the quadratic 
and the quartic parts in terms of the quark cre- 
ation/annihilation operators. For the vacuum energy 
density one has 



— (0\TH[<p]\0) 



d 3 p 



(27T) 



. [A{p) + m ] S(p) 
[B(p)+p}C(p)) : (10) 



where Vol is the three-dimensional volume; the degen- 
eracy factor g counts the number of independent quark 
degrees of freedom, 



g = (2s + l)N c N f , 



(11) 



with s = \ being the quark spin; the number of colours, 
Nc, is put to three, and the number of light flavours, Nf, 
is two. Thus we find that g = 12. 

To arrive at the mass gap equation, we compute the 
Fourier transform of the potential. To regularize the in- 
frared (IR) part of the potential, a modified version of 
the potential ( [3]) [13 is convenient for f3 > — 1, 



V (r) = ~aK fi+1 \rfe- m ^, 



(12) 



where m plays the role of the regulator for the infrared 
behaviour of the interaction, but the limit m — > is 
understood. We get for the Fourier transform, 



y(k) 



a \ d 3 r K 



mr^-ik-r 



- a ^L / drK o 0+1 r 0+1 e- mr sin kr 
k Jo 

(J3 + 2) arctan JM 



2) sin 



-Q 



\k\ 



(k 2 + m 2 )' 



4 7 r J K'o /3+1 r(/3 + 2) sin 



+a 



III 

2 



\k\ 3 +P 



(13) 



where, with the simple density-density harmonic interac- 
tion 15[ , the integral of the potential is a laplacian and 



where the Fourier transform only exists for (3 > — 2. For 
smaller exponents the Fourier transform is UV divergent. 
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For the generalized power-like potential (fl2)) . the angu- 
lar integrals necessary to compute the intermediate func- 
tions A(p) and B(p) are, 



h = 



dio V(k 2 + k' -2u>kk') 



^T(/3 + 2) sin ^ x 



/ +1 
2 

(TT/3) 2W 



\k + k'\ 1+ P \k-k'\ l 



+f3 



h = 



dtu V(k 2 + k' 2 - 2u>kk') w 



-aA-KK Q p+L T{f3 + 2) sin ^ x 



y3+l T 

2 1 

(1 + 0) 2kk' 
4 



1 



1 



\k + k'\ l +' 3 Ik-k'^+P 
1 



-l + /3)(l + /3) {2kk'f 
1 1 



\k + k'\- l +v \k-k'\- l +e 

and we get, 



(14) 



C 

A( P ) 



r(l + /3)sin^ 



B(p) = p + C 




(15) 



where, for the sake of convenience, we continued the in- 
tegral to the negative values of k assuming that m(p) is 
an even function, as it would happen in 1+1 dimensions. 

Consequently, the mass-gap equation (JT)) takes the 
form of a non-linear integral equation for the constituent 
mass m(p) 



m(j>) 



i + / dkC-D(k) 

J —OC P 





pk 


{• 


\p-k\ 1 +f } _ 



1 



\p-k\- l +i 3 



pm(k) — km(p)] 

km{p) > (16) 



III. ANALYTICAL PROPERTIES OF THE 
MASS GAP EQUATION 

We now analyse dimcntionally the mass gap equation, 
and study possible infrared (IR) and ultraviolet (UV) 
divergences. 

In Section II the possible exponents (3 are already lim- 
ited to (3 > — 2, since in eq. (fT3|) the Fourier transform of 
the potential does no exist for smaller exponents. Now 
we address in a dimensional analysis the stability of a 
possible non-trivial vacuum. Let us assume that a solu- 
tion m{k) exists, minimizing the vacuum energy. Then 
we arbitrarily rcscale the solution, 



m(k) — > m(K k) 



(17) 



and, if the vacuum energy does not have an absolute 
minimum for k = 1 , the vacuum is unstable and thus our 
assumption was wrong. We may now simply evaluate the 
vacuum energy as a function of the dimensionless factor 
«. We get, from eq. (TIT 



and this is the main object of our studies. 



£ vac (n) = cift 4 + c 2 m K 3 + c 3 K Q ^ +1 K- p+3 (18) 

the vacuum energy density, with a dimension of the 
fourth power of momentum. In eq. (|18[) . the Cj are con- 
stants, equal to the different integrals in the vacuum en- 
ergy density (| L0[) . Thus, in the chiral limit which is the 
one mattering here, the kinetic energy density scales like 
k 4 . We can show that c\ is positive, and this prevents 
the vacuum to be UV unstable, providing the potential 
term has a smaller scaling power than the kinetic energy 
density. Thus for j3 > — 1 there may be a solution. The 
nicer case is the one of the linear potential where the 
£vo,c(k) has a perfect Mexican hat shape. For (3 < — 1, 
the potential always wins the kinetic term, moreover for 
an attractive potential we can show that the constant C3 
is negative, and thus the vacuum is unstable. For the 
Coulomb case, (3 = — 1 both terms scale equally, and 
there is either a trivial solution m(k) = if the C3 < c\ 
or the vacuum is unstable if c\ > C3. Thus we show that 
there may be a stable and non-trivial m(k) > solu- 
tion to the mass gap equation only for (3 > —1. Since 
the present paper is specialized in negative exponents, 
we are interested in solving the mass gap equation for 
-l</3< 0. 

In the present case of a negative exponent (3, we now 
show that the gap equation (JTJjis IR finite. Notice that 
in the mas gap equation any possible IR divergence may 
only occur in the two denominators with a power of 
\p — k\. The second denominator | p _ fc jL 1+)j is clearly IR 
finite since the exponent 1 — f3 > 1. The first denomi- 
nator | p _ fc 1 |i +/J tends to an IR divergence when j3 — > 0, 
however this divergence is cancelled by the numerator 
pm(k) — km(p). Thus the mass gap equation ([7]) for neg- 
ative exponents is quite different from the mass gap with 
positive exponents, where an exact cancellation of the IR 
divergences of these two different terms would occur, but 
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nevertheless would be technically harder to implement in 
the mass gap equation. 

In the present case of IR finiteness, we now focus in 
the UV sector of the equation. In what concerns the UV 
limit, each separate term in the integrand of the mass 
gap equation ([7]) may be UV divergent in the limit of 
(3 — » — 1, when the integrand momentum k tends to ±oo. 
To study whether the UV divergences cancel, it is conve- 
nient to perform momenta expansions in the integrand. 
There are two different expansions of interest, one where 
the momentum in the integrals is much larger than the 
external momentum, corresponding to the the limit of 
|fc| >> where the mass is not limited. A second pos- 
sible limit is the one of k,p >> m when we are interested 
in large external momenta, and we assume that the mass 
is limited. We start by expanding the integrand in the 

limit when - — > 0. 
p 

In order to utilize as much as possible the cancellations 
of UV divergences of the different terms in the mass gap 
equation, we not only sum all the terms but also return 
to a momentum integral from to oo. Then, if I(k) is 
the integrand of the mass gap eq. (|16|) . for the expansion 



in | of the mass independent terms we get, 

| | [(3 + 0)m(f>) - 3m(*)l + i(3 + /3)(2 + ji) x 



2 4 



[(5 + /3)m(p)-5m(fe)]^ + o(^) 



(19) 



thus, at leading order in ^, the mass gap equation |T 
can be rewritten as, 



aKoW . tt/3 2 

C = — 7) r ( 2 + /3) sm ^-q 



(20) 



m Q = m(p)-J dkC'j^ [(3+j9)m(p)-3m(fc)] ■ 

The eq. (|2"0"|) is UV divergent for (3 < —1 but it is indeed 
UV finite for — 1 < (3 < 0, and thus we can proceed in 
our search of it's solution. 



IV. ALGEBRAIC SOLUTION OF THE MASS 
GAP EQUATION IN THE CHIRAL LIMIT 

We now utilize algebraic methods to solve the mass 
gap equation. We first set the mass gap equation in a 
form close to an eigenvalue equation and show that the 
eigenvalues of this equation are real. The solutions of 
the mass gap equation are the roots of the linearized 
eigenvalue equation, and this provides a fast convergence 
to the solution. We provide the solution in the limit 
where the UV contribution leads the mass gap equation. 



We consider the chiral limit of mo — > 0. We first ad- 
dress the mass gap equation in three momentum dimen- 
sions (3d). The 3d mass gap equation ([7]) van be rewrit- 
ten as 



1 



D{p) 
1 
2 
1 
2 



D(p) m(p) 
d 3 k 



V(p-k)(p-A:) 



kD(k) 



(2tt) 3 



P D(p) 
D{k)m{k) 



D(p) m{p) 



(21) 



and this is similar to a 3d - like algebraic quasi-linear 
equation with a symmetric matrix for a vector D(k)m(k), 
since the integrand only depends on the distance p k. 

To solve the mass gap equation, we start by fixing the 
denominator functions D(k) inside the square brackets 
] of eq. (|2Tj) . with our best initial guess. Then we apply 
the eigenvalue method to the resulting matrix equation. 
Clearly, the eigenvalues are real since the matrix is real 
and symmetric. A solution of the mass gap exists when 
the matrix has a root. In eq. (|2ip we have three matrices. 
The first one is the identity and has positive eigenvalues 
growing with the mass function m(p). The other two 
matrices have eigenvalues with little dependence on the 
mass function m(k). If, in the limit of vanishing m(fc), 
there are one or more negative eigenvalue, then when we 
increase the mass m(k), the first matrix increases and 
eventually it is able to cancel the negative eigenvalues of 
the second plus third matrices. In that case we find the 
desired roots, and we solve the mass gap equation. The 
number of solutions is equal to the number of negative 
eigenvalues of the matrix computed in the masslcss limit. 

To actually solve the mass gap equation it more con- 
venient to solve the radial momentum version (Id) of the 
mass gap equation. For the Id mass gap we can start 
from eq. (|16p . and rewrite it is a symmetric form, 



D(p) 



D(p)pm(p) + 

1 1 

{-l + P) \p-k\- l +P 

-1 



dk^k^ 
P 3 D(p) 

D(p)pm(p) 



pk 



Ip-k^+P 



dkC 



Ip-k^+f 3 



D(k)km{k) 



(22) 



and again particular we have two terms, one fully di- 
agonal, and another explicitely symmetric since it only 
depends on the diagonal distance \p— k\. Thus the matrix 
is hcrmitean, and the eigenvalue equation now applies to 
a vector D(k) km(k). Again, the number of solutions is 
equal to the number of negative eigenvalues of the matrix 
computed in the massless limit of D{k) — ► jr. 

We now solve the mass gap equation in the limit where 
it is lead by the UV contribution. From the UV lead mass 
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FIG. 3: We show, (left) the quark mass m, and (right) log(m), solution of the mass gap equation in the chiral limit, plotted as 
a function of the exponent e, in dimensionless units of Kg = 1. We consider here a = n/12 as in the Luscher term [J] computed 
in static Lattice QCD potentials. 



gap cq. (|20[) we get the matrix equation, 



= 



1 - 



dkC 



D(k) 



dkC' 



D(k) 



(3 + /3) 
m(k) 



(23) 



Let us consider that we discretize the momenta, for in- 
stance in a lattice of N points, where the correct solution 
is found in the limit ./V — > oo. Then the mass gap equa- 
tion (|23]l is a iV x W matrix equation. In this case the 
corresponding matrix is not symmetric, but actually we 
can find the eigenvalues exactly and show that they are 
real, since one of the matrices is a constant and the other 
is a projector. 

In particular, the eigenvalue equation applies to the 
vectors m(k). The first term in cq. |23|) . proportional to 
m(p), is constant, and is thus proportional to the identity 
matrix. The term integrating in m(k) does not depend on 
p and thus it is a projector on a constant vector mi(j>) = 
est. 

This implies that one eigenvector, say Ai(p), of the 
N x N matrix is constant, 



Ai 



dkC ' ( 3 ' 

'Mm* 



0) 



= i- 



\k\ x +p 

aKL |fc|l+/3 P 



(24) 



And, since in the projector all the lines of the discretiz- 
ing matrix are identical, the matrix has N — 1 linear de- 
pendences, and this implies that all the other TV— 1 eigen- 
vectors Xi(p), for i > 1, are cancelled by the projector, 
and thus all their eigenvalues are identical, simply given 



by the matrix proportional to the identity. The other 
eigenvalues Xi(p) , i > 1 cancelled by the projector are, 



dfcC'jjpL(3 + /3)} . (25) 



Now, notice that C < if we consider an attractive, 
i.e. negative shorter range potential with —a < 0, and 
with a negative (3 ~ — 1, and leading us to, 



Ai-1 = - 



A* - 1 



dkC 



k^+p 

\k\ 



1+/3 



(3 + 0) > 1 . (26) 



Thus we can have one, and only one root. The other 
eigenvalues are positive (and quite large if the integral is 
nearly UV divergent). Thus this differs from the confin- 
ing potentials where a whole tower of replicas was found 
by Bicudo and Nefediev 

We now concentrate on the root to determine the 
constant eigenvector m that solves the mass gap equa- 
tion (|20[) . It is convenient to rename the exponent to 
P = — 1 + e, where we are interested in the exponent 
range of < e < 1. The mass gap equation is then, 



Ai 



1 - 



(27) 



dkC 



I -C {-! + €)- 



Vk 2 +m 2 \k\' 

e rmr(f) 

20F 
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and the solution is, 



r(^)r(§; 



2^ 



= Kn 



r(l + e) 7r(l-e)l-er 
a — sin ■ 



2- 



2 f 



a ^2 + 37 + ^(5) 
•ire 



()7T 

2 + 3, + ^) 

» ' +0 



o(e) 
0) 



(28) 



and this diverges very fast in the limit of vanishing e, but 
for finite e it occurs that the negative (3 exponents indeed 
produce chiral symmetry breaking. 

The solution of eq. ([28]) for the dynamical quark mass 
is plotted in Fig. [3] as a function of e. Wc find that a 
small but finite mass m ~ 0.05Kq, almost independent of 
e. in the range 0.15 < e < 0.9. For e ~ 1, corresponding 
to (3 ~ the mass vanishes. The mass explodes for e ~ 0, 
close to the coulomb potential. 



THE QUARK ENERGY, IN THE CHIRAL 
LIMIT AND WITH A BARE MASS 



To compute the quark energy, it is first convenient to 
write the mass gap equation |(7J) as, 



Mp) = b(p) 

m(p) p 



and then we get for the quark energy 



(29) 



E(p) 



D{p) [m(p)A(p)+ P B( P )} 
A(p) 



y / m{p) 2 + p' 2 



m(p) 



\/ m 2 + 



1 - e 



(30) 



where to is computed in cq. (J28J). 

Since to is UV divergent when e — > 0, In order to have 
a limited quark mass when e is small, one needs to in- 
clude in the mass gap equation, as a counter term, a 
bare quark mass mo- I n this case we may use a different 
limit, where to «,p, k, still considering for simplicity a 
constant m(p) = m, and we get for the A(p) and B(p) 



integrals, 
A(p) 



rn 



C / k 2 dk 



m 

toq C 2 



{2pk 



Ip-k^+P 



to p e 1 — e 



B(p) 



m 



1 + C 



1 1-7-loghfc 



o(e) 



2pfc 



1 



1 



(-1 + /3) (2pfc) 2 
C 4 



1 



2a 

1 + T 



p- (l- e )(3-e) 

J_ 4-3 7 -31og( ^ 
7re 



(31) 



3tt 



o(e) 



Thus, assuming an approximately constant to we get for 
the bare mass too which acts here as rcnormalization 
counter term, 



to 



a 
3 7r e 



0(e°) 



(32) 



and thus, if we specialize in the minimal subtraction 
scheme, wc get a finite mass. Nevertheless the quark 
energy remains UV divergent in the limit of e — > 0, 



E(p) = y/P 2 



^— +0 e 1 

37T6 



(33) 



A further rcnormalization of the quark energy may be 
necessary, but only upon studying the Bethe Salpeter 
equation for mesons, which goes beyond the scope of 
the present paper. Different methods to renormalize the 
Coulomb potential in the Coulomb gauge have been ap- 
plied by the Zagreb group [13] and by Szczepaniak and 
Swanson [27[ utilizing the Glazek and Wilson method 

M. 



VI. CONCLUSION 

We study xSB driven by power-law potentials with 
a negative exponent f3 < 0. These potentials are non- 
confining. We extend a previous study performed for 
confining potentials with a positive exponent [l| . In that 
study, chiral symmetry already vanishes when f3 — » 0, 
nevertheless, since the Coulomb potential is negative, it 
is natural to reverse the sign of the potential, and then 
%SB may occur for f3 < 0. 

We work in momentum space, and the existence of a 
Fourier transform is limited to p > —2. Then, with a 
dimensional analysis, we show that we may only have 
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a stable \ S B vacuum for /? > — 1. Thus we study in 
detail the non-confining potentials with — 1 < j3 < 0, i. e. 
in the exponent range limited by the Coulomb potential 
and the logarithmic potential. 

We find that %SB also occurs for the studied negative 
exponent power-law potentials generating dynamically a 
finite quark mass to. 

Moreover we find qualitative differences to the mass 
generated with confining potentials. First, we find one 
and only one non-trivial solution of the mass gap equa- 
tion, whereas for confining potentials an infinite tower 
[l[ of false, excited vacua, sometimes called replicas, is 
found. Also, the solution found has an approximately 
constant mass, as in Fig. [3j i. e we find no IR enhance- 
ment of the quark mass, whereas the confining potentials 
studied previoulsy [l[ produce a significant IR enhance- 
ment of the quark mass, as in Fig. [2] 

We also find an UV divergence of the dynamical mass 



to, in the Coulomb potential limit of the exponent (3 — > 
— 1. This is consistent with the well known necessity 
to apply a renormalization program when the Coulomb 
potential is used. 

This work may be a starting point for the study of xSB 
at the deconfincmcnt transition, say with finite T, finite 
fj, or large Np. 
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